Parton rescattering effect on the charged hadron forward-backward multiphcity 

correlation in pp colhsions at y^=200 GeV 
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The parton rescattering effect on the charged hadron forward-backward muhipUcity correlation 
in pp collisions at ^/s=200 GeV is studied by a parton and hadron cascade model, PACIAE, based 
on the PYTHIA model. The calculated multiplicity and pseudorapidity distribution of the final 
state charged hadron are well compared with experimental data. It turned out that the final state 
charged hadron pseudorapidity distribution are different from the initial state charged partons. The 
parton rescattering effect on the charged hadron forward-backward multiplicity correlation increases 
with increasing parton rescattering strength in the center pseudorapidity region (I77I < 1). However, 
this effect becomes weaker in the outer pseudorapidity region {\r]\ > 1). 

PACS numbers: 24.10.Lx, 24.60.Ky, 25.75.Gz 



I. INTRODUCTION 

The study of fluctuations and correlations has been 
suggested as a useful means to reveal the mechanism of 
particle production and the formation of Quark-Gluon- 
Plasma (QGP) in the relativistic heavy-ion collisions 
[U, 0]. Correlations and fluctuations of the thermody- 
namic quantities and/or the produced particle distribu- 
tions could be significantly altered when the system un- 
dergoes phase transition because the degrees of freedom 
is largely different between the hadronic matter and the 
quark-gluon matter. 

The experimental study of fluctuations and correla- 
tions becomes a hot topic in relativistic heavy ion colli- 
sions with the availability of high multiplicity event-by- 
event measurements at the CERN-SPS and BNL-RHIC 
experiments. An abundant experimental data have been 
reported [1, i, H, i, 0, i, i, M, [U where a lot of new 
physics are arisen and urgent to be studied. 

A lot of theoretical investigations have also been re- 
ported such as Kefs. [H [ll [H, [H, [3. We have 
used the PYTHIA model to investigate the strength of 
charged hadron forward-backward multiplicity correla- 
tion in pp and pp collisions at y^=200 GeV [16]. It ar- 
gued that a factor of 3-4 apparent discrepancy between 
UA5 pp data [13] and STAR pp data [11] can be at- 
tributed to the differences in detector acceptances and 
observed bin interval in both experiments. In [l5| the 
back-to-back parton scattering was considered as the ori- 
gin of final state hadronic covariance. They assumed 
the details of hadronization are not essential and related 
the back-to-back parton scattering angles to a Gaussian- 
like hadronization function. Then they derived the fi- 
nal state charged hadron forward-backward multiplicity 



covariance. Without more dynamical inputs^their re- 
sults were well compared with STAR data [ig], "thus 
dispelling the notion that correlation length has any fun- 
damental significance" . Stimulated by this interesting 
conclusion, a parton and hadron transport model, PA- 
CIAE [1^, is employed in this paper to investigate the 
effect of parton rescattering in parton evolution stage on 
the final state charged hadron forward-backward multi- 
plicity correlation in the pp collisions at \/s= 200 GeV. 

Following [23] the strength of charged particle forward- 
backward multiplicity correlation, b, is defined as 



(n/nb) - (n/)(nb) 



(1) 
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where Uf and Ub are, respectively, the number of charged 
particles in forward and backward pseudorapidity bins 
defined relatively and symmetrically to a given pseudo- 
rapidity rj. The (ny) ((jib)) is the mean value of (rib), 
the Z3jj, and Dj:j are the forward-backward multiplic- 
ity covariance and forward multiplicity variance, respec- 
tively. One always studies 6 as a function of the center 
distance of two forward and backward pseudorapidity 
bins (Ary) and the rj acceptance. 



II. THE PACIAE MODEL 

The parton and hadron cascade model, PACIAE [l9| . 
is based on PYTHIA Jl| which is a model for high en- 
ergy hadron-hadron (hh) collisions. The PACIAE model 
is composed of four stages: parton initialization, parton 
evolution (rescattering) , hadronization, and hadron evo- 
lution (rescattering). 

1. PARTON INITIALIZATION: In the PACIAE 
model, a nucleus- nucleus collision is decomposed into 
the nucleon-nucleon (NN) collisions based on the col- 
lision geometry. The NN collision is described with the 
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PYTHIA model, i.e. it is decomposed into the parton- 
parton collisions. The hard parton-parton collision is de- 
scribed by the lowest-leading-order (LO) pQCD parton- 
parton cross section [2^ with modification of parton dis- 
tribution function in the nucleon. And the soft parton- 
parton interaction is considered empirically. The semi- 
hard, between hard and soft, QCD 2^2 processes are 
also involved in the PYTHIA (PACIAE) model. Be- 
cause the initial- and final-state QCD radiation added to 
the parton-parton collision process, the PYTHIA (PA- 
CIAE) model generates a partonic multijet event com- 
posed of quark pairs, diquark pairs and gluons for a 
NN [hh) collision. That is followed, in the PYTHIA 
model, by the string-based fragmentation scheme (Lund 
string model and/or Independent Fragmentation model). 
Thus a hadronic final state is reached for a NN {hh) 
collision. However, in the PACIAE model above frag- 
mentation is switched off temporarily, so the result is 
a partonic multijet event instead of a hadronic state. If 
the diquarks (anti-diquarks) are split forcibly into quarks 
(anti-quarks) randomly, the consequence of a NN (hh) 
collision, thus a nucleus- nucleus collision, is its initial 
partonic state composed of quarks, anti-quarks, and glu- 
ons. 

2. PARTON EVOLUTION: The next stage in the 
PACIAE model is the parton evolution (parton rescat- 
tering). Here the 2^2 LO-pQCD differential cross 
sections [l^ are employed. The differential cross section 
of a subprocess ij — > kl reads 



3. HADRONIZATION: The parton evolution stage 
is followed by the hadronization at the moment of par- 
tonic freeze-out (no more parton collision at all). In 
the PACIAE model, the phcnomcnological fragmenta- 
tion model and coalescence model are supplied for the 
hadronization of partons after rescattering. The String 
Fragmentation (SF) model is adopted in this paper. We 
refer to [11] for the details of the hadronization stage. 

4. HADRON EVOLUTION: After hadronization the 
rescattering among produced hadrons is dealt with the 
usual two-body collision model. We neglect the hadronic 
rescattering in pp collisions as usual. The details of 
hadronic rescattering see (23j . 



TABLE I: Charged particle multiplicity {\r)\ < 5.4) in pp 
collisions at y^=200 GeV. 





PACIAE 


Exp. data^'^ 




parton 


hadron 


no parton scat. 


5.96 


21.2 (17.7)* 




18% parton scat. 


6.13 


21.7 (18.5) 


19.9 ± 2.2 


58% parton scat. 


6.15 


21.9 (18.7) 





^ after parton rescattering, in full 77 phase space. 

^ in final hadronic state. 

^ inelastic pp collision data taken from [24]. 

* value given in bracket is calculated for inelastic pp collisions. 



7 — = i\ 

dt s 



(2) 



where the K factor is introduced for considering the 
higher order pQCD and non-perturbative QCD correc- 
tions as usual and as stands for the effective strong cou- 
pling constant. Taking the process qiq2 qiq2 as an 
example one has 
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where the s, t, and u are the Mandelstam variables. 
Since it diverges at t=0, it has to be regularized by in- 
troducing the parton colour screen mass fj, as follows 



E 



9{t- fj.^y 



(4) 



9192— » 91 92 

The total cross section of the parton collision i+j then 
reads 



)-E 



k.l 



dt- 



ij — ^kl 



dt 



(5) 



With above total and differential cross sections the par- 
ton evolution (parton rescattering) can be simulated by 
the Monte Carlo method. 



III. PARTON RESCATTERING EFFECT ON 
HADRON MULTIPLICITY CORRELATION 

As we aim at the physics behind the experimental data 
rather than reproducing the data, model parameters are 
fixed in the calculations. All calculations are for y^=200 
GeV Non-Single Diffractive (NSD) collisions except that 
marked especially. In order to see the parton rescattering 
effect in the parton evolution stage, we design no, weak, 
and strong parton rescattering cases as follows: 

• No parton rescattering: K=0, without parton 
rescattering at all. 

• Weak parton rescattering: K=l and 11^— OA 
GeV^/c^, with nearly 18% charged partons par- 
ticipating the rescattering. 

• Strong parton rescattering: K=3 and im^=0.1 
GeV^/c'*, with nearly 58% charged partons par- 
ticipating the rescattering. 

Table |T] shows that although we did not adjust the 
model parameters the PHOBOS multiplicity data of fi- 
nal state charged hadrons 2J| are well reproduced. From 
the PACIAE results in this table one sees that the multi- 
plicity of initial state charged partons {u+d+s and their 
anti-particles) and final state charged hadrons increase 
weakly from no, to weak, and to strong parton rescat- 
tering case. This increase in percentage is not as strong 
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FIG. 1: (a) The pseudorapidity distributions of the initial 
state charged partons and final state charged hadrons and 
(b) the parton rescattering effect on the final state charged 
hadron pseudorapidity distribution in pp collisions at y'^— 
200 GeV. The PHOBOS inelastic ri distribution data are 
taken from [23. 



as the increase in percentage of participating charged 
parton given in the rescattering case definition. This is 
because the total number of initial state charged partons 
is just around 6 and only the 2 —^ 2 parton rescattering 
processes are considered in the parton evolution stage. 

The calculated pseudorapidity distributions of initial 
state charged partons and final state charged hadrons 
are compared in Fig. [TJa). The full and open squares 
are, respectively, the PHOBOS inelastic pp collision data 
(taken from |24| ) and the corresponding theoretical re- 
sults for final state charged hadrons. And the open 
circles are theoretical results calculated for final state 
charged hadrons in NSD pp collisions. One sees that the 
PHOBOS data agree well with the NSD calculations but 
are higher than the inelastic calculations in the center rj 
region. This is consistent with the fact that the charged 
hadron multiplicity in inelastic calculations is lower than 
the PHOBOS data as shown in Tab. [H Comparing the 
rj distribution of final state charged hadrons to the one 
of initial state charged partons (open triangles) we know 
that the parton rescattering and fragmentation fill up the 
wide valley between two fragmentation peaks at I77I ~ 5 
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FIG. 2: The parton rescattering effect on the final state 
charged hadron forward-backward multiplicity correlation 
strength in (a) central pseudorapidity range (|?7| < 1) and 
(b) outer pseudorapidity range (1 < \ri\ < 4) in pp collisions 
at -^5=200 GeV. The open squares, open circles and open 
triangles are for case of no, 18% and 58% participant parton, 
respectively. The experimental data are taken from p^ . 



in the T] distribution of initial state charged partons. 

Fig. [ijb) shows the parton rescattering effect on 
the pseudorapidity distribution of final state charged 
hadrons. One sees that in the I??! < 4 region the pseudo- 
rapidity distributions move upward monotonously from 
no, to weak, and to strong parton rescattering case. 
However, in the outer region one sees nearly the opposed 
situation. 

Figure [2] shows the parton rescattering effect on the 
forward-backward multiplicity correlation strength of fi- 
nal state charged hadron in pp collisions at y/s=200 
GeV. The corresponding experimental data [l8| are given 
by full squares. We see in panel (a) that the correla- 
tion strength increases with increasing strength of par- 
ton rescattering in the central pseudorapidity region of 
I77I < 1 (or A77 < 2). While panel (b) shows that the par- 
ton rescattering effect becomes weaker in the outer pseu- 
dorapidity region \r]\ > 1 (or Ar] > 2). It has to point 
out here that the rj acceptance in panel (b) is nearly four 
times larger than the one in panel (a) so the b at A77 = 2 
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FIG. 3: The parton rescattering effect on the statistical and 
dynamical correlations of final state charged hadron in central 
pseudorapidity range < 1) in pp collisions at ^/s—200 
GeV. 



in panel (b), for instance, is larger than the correspond- 
ing one in panel (a). 

We also use the mixed events method [16;] to study 
the parton rescattering effect on the statistical and dy- 
namical (non-statistical) correlations separately. The re- 
sults are shown in Fig. [3] We see that the initial state 
parton rescattering have almost no effect on statistical 
correlations. That is because the statistical correlation 
is steaming from the multiplicity fluctuation in the de- 
tected 77 bins and the parton rescattering have only little 
effect on the multiplicity of final state hadrons as shown 
in Tab. [D However, one sees that the dynamical corre- 
lation strength is increasing monotonously from no, to 
weak, and to strong parton rescattering case in central 
pseudorapidity range (I77I < 1). That shows the parton 
rescattering increase the dynamical correlation among 



partons, it survived the hadronization, and effected on 
the dynamical correlation of final state hadrons. 



IV. SUMMARY 

Using a parton and hadron cascade model PACIAE 
we have investigated the parton rescattering effect on 
forward-backward multiplicity correlation strength of 
the final state charged hadron in pp collisions at ^/s = 
200 GeV. The calculated multiplicity and pseudorapidity 
distributions of the final state charged hadron are well 
compared with the corresponding experimental data, ft 
turned out that the pseudorapidity distribution of final 
state charged hadron is different from the initial state 
charged parton. The parton rescattering effect on the 
forward-backward multiplicity correlation of the final 
state charged hadron increases with increasing parton 
rescattering strength in the center pseudorapidity region 
(1 77! < 1). This increase becomes weaker in the outer 
pseudorapidity region {\r]\ > 1). However, the final state 
hadron correlation strength is related to the initial state 
partons. 

The parton rescattering effect on the forward- 
backward multiplicity correlation of final state charged 
hadron is not as strong as expected from the increasing 
of participating charged parton. That is because the av- 
erage total number of initial state charged partons is just 
around 6 in above collisions and the parton rescatterings 
considered in parton evolution stage are all 2 — > 2 pro- 
cesses. This is consistent with the notion that the small 
probability of QGP formation in pp collisions at RHIC. 
One may expect much stronger parton rescattering ef- 
fect in the nucleus- nucleus collisions at the same energy 
and even in the pp collisions at LHC energy. 
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